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Abstract 

The differences between two standard Monte Carlo models, Lepto and Ari- 
adne, for deep inelastic scattering at small-x is analysed in detail. It is shown 
that the difference arises from a 'unorthodox' suppression factor used in Ari- 
adne which replaces the normal ratio of parton densities. This gives rise to 
a factor that qualitatively is similar to what one would expect from BFKL 
dynamics for some observables like the energy flow and forward jets but not 
for the 2+1 jet cross-section. It is also discussed how one could use the 2+1 
jet cross-section as a probe for BFKL dynamics. 
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Deep inelastic electron proton scattering at small-x (with x being the ordinary Bjorken- 
x) has drawn much theoretical and experimental interest the last couple of years. Since the 
observation [|I],0 of the predicted rise of the deep inelastic structure function F2 at small- 
X, the main question has been whether the 'unconventional' so called BFKL dynamics 
can be observed or whether the interactions can be described by ordinary leading log(Q^) 
dynamics in perturbative QCD as given by the DGLAP |^ evolution equation. From the 
inclusive measurement of F2 it seems not possible to draw any such conclusions based on 
the presently available data and it is also theoretically questionable whether it will be at all 
possible as long as one only considers one observable 

Other observables that have been proposed to see effects of BFKL dynamics are the 
transverse energy flow and forward jets [If. One problem in looking for possible effects of 
BFKL dynamics is that so far there has only been analytical calculations on parton level and 
the hadronisation effects are difficult to estimate. In this respect the forward jets analysis 
is more promising since the hadronisation corrections are expected to be small for a jet 
cross-section. However, the forward jets are close to the target remnant jet which makes it 
difficult to disentangle the two. 

So far there exists no Monte Carlo (MC) model based on the BFKL equation even though 
there has recently been significant progress in creating such a model ^ based on the CCFM 
equations [1^] which interpolate smoothly between the BFKL and DGLAP equations. In the 



mean time the Ariadne |jTT]] MC which is based on the colour dipole model for DIS [|T2| has 
often been used to estimate the possible effects of BFKL dynamics. The argument for this has 
been that it also contains parton emissions which are non-ordered in transverse momentum. 
However, the present paper will show that the difference between Ariadne and ordinary 
DGLAP evolution is mainly due to a 'unorthodox' suppression factor used for additional 
emissions from dipoles connected with the proton remnant. Comparing the results from the 
Ariadne MC with data has also shown rather good agreement for the transverse energy 
flows and forward jets whereas a pure DGLAP Monte Carlo fails. There are, however. 



large hadronisation uncertainties and it has been shown that the Lepto |T3| MC model 
based on leading order (LO) QCD matrix-elements and leading logarithmic corrections in 
the form of parton showers according to the DGLAP evolution equation supplemented with 
non-perturbative hadronisation effects can also describe the data on transverse energy flow 

a- 



An observable which at first sight might seem more or less uncorrelated to the question 
of whether one observes BFKL dynamics or not is the 2+1 jet cross-section. Since this is 
a genuinely hard process it should be well described by conventional perturbative QCD as 
long as the more complicated forward region is excluded. However this is not the case. As 
argued in this paper the 2+1 jet cross-section should exhibit features of BFKL dynamics 
when the mass s of the jet-system becomes much larger than the momentum transfer t 
from the incoming parton. In other words one should be able to observe effects of BFKL 
dynamics as the rapidity difference {Ay = ln{s/\t\)) between the two jets increases. This is 
just the same kind of mechanism as in dijet production in hadron collisions as proposed by 
Mueller and Navelet |15|. There has also been a proposal [jl6| to use the 2+1 jet cross-section 



for small jet-systems as a probe of BFKL dynamics since this would essentially probe the 
x-distribution of gluons at small x where effects of BFKL dynamics should be possible to 
observe. A complication with this proposal is that the gluon density is not an observable 
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and the constraints from data on what it should be in a DGLAP evolution scenario are small 
especially if one takes into the account the theoretical uncertainties in defining the gluon 
density. There are also experimental difficulties in finding jets for small mass jet-systems. 

At large values of photon virtuallity (typically > 100 GeV^) the measured jet 
cross-section has been shown to be in good agreement p]7|-[20[| with the Monte Carlo model 
Lepto and next-to-leading order (NLO) perturbative QCD calculations using programs 
such as DisJET and Pro JET [^]. In addition, the Ariadne MC also describes the 
Vcut dependence (where i/cut is the cut used in the JADE jet- definition, mjj < UcutW^) of 
the data with a similar accuracy to the one found for Lepto whereas the Q^-dependence is 



not so well described |jT8|,IT9[. These cross-section measurements for large Q have also been 



used for the extraction of as and the observation of its running |ll9| , p0 |. 



However, if one instead looks at the small-x and small region, where x = Q^/2P ■ q 
with q and P being the photon and proton four momenta respectively, but still avoiding the 
more complicated forward region, the picture does not seem to be so clear anymore. There 



is still fairly good agreement between Lepto and data [18,19| but the two MC models give 



quite different results. As an example Table | gives the 2+1 jet cross-section on matrix- 



element level as calculated with Lepto 6.5, Ariadne 4.08 and Mepjet 1.1 |^ which is 
a NLO perturbative QCD calculation. The kinematical region used is: 0.0001 < x < 0.001, 
0.05 < Hb < 0.7 [i/B = q ■ P/k ■ P with k being the momentum of the incoming lepton), 



5.0 < < 70.0 GeV^ and the MRS A parton distributions pil have been used. The jets 



have been defined using the cone- algorithm in the hadronic center of mass system (hems) 
with a cone-size of i? = ■\/Arp~+~A0^ = 1 and requiring the jets to have transverse momenta 
p± > 5 GeV and 0.1 < z < 0.9 (where z = P-ji/P-q with jj being the jet-momenta) to avoid 
the forward region. This gives jets which have pseudo-rapidities in the photon-hemisphere 
in the hems and —2 < 77 < 2 in lab, thus avoiding the forward region. The cut in z also 
decreases the effects from parton showers on the 2+1 jet cross-section to less than 10 percent 



TABLE I. Jet cross-sections for 0.0001 <x < 0.001, 0.05 < vb < 0-7 and 5.0 < < 70.0 GeV^ 

Program ;u|., /i^ a [nb] 

Ariadne Q^Px 6^38 

Lepto 152 

Mepjet LO 1.49 

Mepjet LO Q'^,{T,p±/2)'^ 1.19 

Mepjet LO {'^P±/2f 1.45 

Mepjet NLO 0.89+0.06 

Mepjet NLO (Sp_l/2)2 0.95+0.04 



The first thing to note is that the 2+1 jet cross-section predicted by the Ariadne MC 
is about four times larger than the one from Lepto. Comparing with the LO result from 
Mepjet using as renormalisation and factorisation scale shows good agreement with 
Lepto whereas the LO result from Mepjet with the same factorisation and renormalisation 
scales as in Ariadne differs a factor five from the Ariadne result. Clearly the LO 2+1 jet 
cross-section in Ariadne does not agree with the result from LO perturbative QCD. The 
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second thing to note is that the NLO corrections are quite large and that they in fact are 
negative such that the NLO cross-section is smaller than the LO one. Of course the question 
of choosing an appropriate renormalisation scale and factorisation scheme and scale has to 
be addressed if one wants to compare the NLO cross-sections with data but for the present 
purposes, i.e. to see that the NLO-corrections cannot increase the cross-section by a large 
factor it is not necessary to pursue these theoretical uncertainties further (for a discussion 
of the renormalisation scale uncertainty see for example ||25||). 




FIG. 1. The Zg, Xp and p±_ distributions from the LO matrix-elements in Ariadne 4.08 (full 
line), Lepto 6.5 (dashed line) and the modified version of Lepto (dotted line) for the two basic 
processes QCD-compton (QCDC) and Boson Gluon Fusion (BGF). 



To see in more detail what the difference between Ariadne and LO perturbative QCD 
is (as represented by Lepto), it is instructive to look at the z, Xp and p_L distributions 
which are given in Fig. |l]. Here Xp = x/E, where ^ is the longitudinal momentum share of 
the parton entering the hard interaction. For clarity the cross-section has been divided into 
two different parts, QCD-compton (QCDC) and Boson Gluon Fusion (BGF), depending on 
the underlying hard process as illustrated in Fig. 0. As can be seen from Fig. |I|, the BGF 
part agrees well between the two programs taking into account the different renormalisation 
scales, but for the QCDC part the difference is dramatic and requires some other explanation. 

In fact, this difference should come as no surprise. The colour dipole model for DIS does 
not use the pure LO perturbative QCD matrix-element for the QCDC process but rather 
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FIG. 2. Illustration of the two leading order 2+1 jet processes, QCD-compton (a) and Boson 
Gluon Fusion (b). 



an approximation. For a given x and the cross-section can (assuming single photon 
exchange) be written as, 

dcrQCDcix,Q'^,z,Xp) _ Cpasfqi^^Q^) , ^2 \'^^ (-w 

where fq is the sum of the quark densities in the proton and g{x,Q'^, z,Xp) is a simple 
function which is not needed for the present discussion. In the colour dipole model for DIS 
the ratio of the parton densities is replaced by a suppression factor due to the extended 
proton remnant, 

where y is the rapidity of the emitted parton in the hems and W is the mass of the hadronic 
system {W^ = (P + qY)- This limits the transverse momentum p± of the emitted gluon 
which corresponds to that only a fraction {fi/p±)°' of the proton remnant takes part in the 
emission. Normally one uses a = 1 but other choices are also possible. Now, for small-x the 
suppression factor used in the colour dipole model for DIS starts to differ substantially from 
the ratio of the parton densities. In the region of interest x is very small and ^ is moderately 
small (typically 10~^ — 10~^) so that this ratio goes essentially as Xp'^^ if one assumes that 
fq{x) oc x~^~'^. For recent parton distributions A is in the order of 0.2 to 0.3 so this gives a 
strong suppression of small Xp in the LO perturbative QCD formula. 

To see that it really is this suppression factor that gives the difference a "toy- model" ver- 
sion of Lepto has been constructed. In this modified version of Lepto the QCDC matrix- 
element has been multiplied by a phenomenological factor x"'' where b = d — O.OSlnl/xp 
with d being a function of x. The value of d was obtained by a fit to the ratio of daqcDC I dxp 
from Ariadne and Lepto for fixed x and Q^. The resulting powers turn out to be easily 
parameterised hy d = 0.2(1 — Xog^^x) in the region of interest (10~^ < x < 10~^). As 
can be seen in Fig. |T]this "toy- model" reproduces the distributions from Ariadne to a very 
good approximation, especially the p^-distribution. Thus, the difference between the QCDC 
cross-section in Lepto and Ariadne can be explained by a factor x"^ which results from 

using the suppression factor 6 [ ^yj^i^pj^Ye-y - P±) instead of fq{^,Q'^)/fq{x,Q'^). As an 
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aside, it is interesting to note that for the BGF-part of the jet cross-section, the suppression 
factor used in Ariadne is indeed the corresponding ratio of the parton densities. 

The important thing to note is that all emissions in Ariadne, where the proton remnant 
is part of the dipole, are treated in the same way as the first emission in the QCDC matrix- 
element. To be more precise the probability for an extra emission, when the proton remnant 
is part of the dipole, is given by 

rfPoc^ef pY-^dy. (3) 

This means that the ratio of parton densities fa{xa,Q'^)/fb{xb,Q^) does not enter as they 
would in a traditional backwards evolution initial state parton shower where, 

4:as fa{xa,Q'^)dpl 

Instead the boundary condition present from the proton is taken into account by the extra 
cut-off in transverse momenta for emitted gluons due to the extendedness of the proton 
remnant. A direct comparison with a traditional initial state parton shower is complicated 
by the fact that the emissions in the colour dipole model for DIS are not ordered in x. But 
even so, it is evident that the difference in the suppression factor gives a large effect for 
the rest of the emissions just as for the QCDC matrix-element. There are two main effects. 
First of all the probability for an emission is increased and secondly the emissions become 
harder, in the sense that they get higher p±, just in the same way as the QCDC part of the 
2-1-1 jet cross-section increases and the pj_ spectrum becomes harder (see Fig. I^c). 

To estimate the magnitude of the effects on the transverse energy flow due to the 'un- 
orthodox' suppression factor the modified version of Lepto with the increased QCDC cross- 
section has been used to calculate the transverse energy flow. The results are shown in Fig. |^ 
where data from HI are compared with Ariadne, Lepto and the modified version of 
Lepto. One should here note that to be able to compare with more or less pure DGLAP 
dynamics the soft colour interactions and the new seaquark treatment in Lepto have been 



shut off for both Lepto versions (see |]T3 for more details) 



As can be seen from Fig. ^ the modified version of Lepto interpolates nicely between 
the Ariadne result in the photon hemisphere where the first emission is important and the 
DGLAP-Lepto result in the proton hemisphere where higher order emissions are important. 
In other words, the modified matrix-element increases the transverse energy flow to the same 
level as in Ariadne in the photon hemisphere, whereas in the proton hemisphere there is no 
difference between the two versions of Lepto since the same initial state parton shower is 
used. It would also be possible to modify the initial state parton shower in Lepto in the 
same way as the matrix-element (essentially by multiplying the splitting function Pq^qdz) 
with a factor z~^) which would lead to an increased transverse energy flow in the proton 
hemisphere. However, already from the result based on the modified QCDC matrix-element 
one can conclude that the difference in transverse energy flow between DGLAP-Lepto and 
Ariadne is mainly due to the 'unorthodox' suppression factor used in Ariadne instead of 
the ratio of parton densities and not due to the difference in p^-ordering. 

The following question then arises. Has this difference in suppression factor anything 
to do with BFKL dynamics? As it turns out, the answer depends on which observable one 
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is interested in. Starting with the 2+1 jet cross-section one first notes that the relevant 
scales for the onset of BFKL dynamics are the momentum transfer t from the incoming 
gluon and the mass s of the jet-system. One expects effects of BFKL dynamics when the 
rapidity difference Ay = ln(s/|t|) becomes large so that a^Ay ~ 0{1) which indicates that 
resummation is necessary. In the hems the relevant ratio becomes \t\/s ^ p^/s = z{l — z). 
This is quite different from the increase with decreasing Xp = Q"^/ (Q'^+s) given by Ariadne. 
One also notes that typical values are of the order of 10~^ — 10~^ which is not very small and 
therefore one does not expect any effects from the BFKL dynamics in the parton densities. 
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FIG. 3. The transverse energy flows for different x and values in the hadronic cms comparing 
Ariadne (full line), Lepto with pure DGLAP (dashed line) and the modified version of Lepto 
(dotted line) with data from HI [^. (The plot has been generated using the HzTool package |p8|). 



The dominating diagram for 2+1 jet production with enhancement from the resummation 
of soft gluons according to the BFKL prescription is depicted in Fig. ^ where the two partons 
from the quark box form one jet. In addition to this process there is also a similar one 
where the incoming and outgoing gluon is replaced by a quark. At sufficiently large rapidity 
difference Ay the contribution from this kind of diagrams should be dominating thanks to 
the t-channel gluon exchange. This is the same kind of diagram that has been analysed for 
the forward jet production ^ but with the important difference that two jets are required. 
So instead of having k± of the forward jet being of the same order as Q for the BFKL 
dynamics to be valid, two jets with the same p± are required. In a sense this is closer to the 



original proposal by Mueller and Navelet |15| to look for BFKL effects in dijet-production 
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FIG. 4. Diagrams with enhancement from BFKL dynamics in the large Ay hmit: (a) lowest 
order diagram (b) all orders diagram with the reggeized gluon indicated with a thick line. In 
addition there is a similar diagram where the incoming and outgoing gluon is replaced by a quark. 



in hadron collisions. 

In a BFKL calculation of the 2+1 jet cross-section one would take the large Ay = ln(s/|t|) 
limit of the 2+1 jet cross-section and reggeize the t-channel gluon which would give a cross- 



section of the type |29 



da 



d^dAy 



oc K 



da 



BFKL 



d^dAy 



largeAy 



(5) 



where the effective parton density PD[, fg{^,iJ?) + ^fq{^,jj?), reflects the domination of 



t-channel gluon exchange for large Ay. The factorisation scale is arbitrary but should be 
of the order of transverse momentum cut-off for the jets, p^,min- 

In a leading order calculation the enhancement factor is equal to one so by forming the 
ratio between the BFKL calculation and the fixed order calculation for fixed ^ one isolates 
the KBFKL-^^ctoT. Thus, as a probe for BFKL dynamics, it should be possible to use the 
2+1 jet cross-section as a function of the rapidity difference Ay between the jets in the hems 
(for fixed ^) and see whether the data start to deviate exponentially from the NLO fixed 
order perturbative QCD predictions when Ay becomes large. To be able to do a quantitative 
analysis one would of course need to actually calculate the proposed cross-section. 

Integrating over t from — ^^j^^ the asymptotic form of the enhancement factor will be 



2| 



K 



BFKL 



exp(asC^Ay4 ln2/7r) 
7asCACzAy/2 



(6) 



Using the standard value for a;sCA41n2/7r = 0.5 this gives. 



Kbfkl oc exp(a^CAAy41n2/7r) 



0.5 



,mm / 



(7) 
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For forward jets —i ~ Q'^ is required and thus Ay ~ ln(l/xp) such that Kbfkl oc (^-^ 



, 0.5 



which also is consistent with the results from 31 



Returning to the cross-section for the QCDC process used in Ariadne, it can be written 
in the following way, 

a,„,(x,Q2) 2rr f,{x,Q^) A|g!) 9[x,L^ ,z,Xp) ^^az 

- ~^ — // n2\ — Q , z, Xp)^dz (8) 

27r fg{x,Q^) \Xp) Xp 

with h being in the order of 0.5 — 1.0 for small x. Qualitatively this gives a ii"-factor x~^ 
which is similar to the one in the forward jet BFKL calculation [^. This also explains 
why Ariadne agrees quite well with data on forward jets just in the same way as the 
forward jet BFKL calculation |^ agrees quite well with data. 

Even though there is some qualitative agreement between the colour dipole model and a 
BFKL calculation for forward jets there are quantitative differences which are important: (i) 
the colour dipole model has an enhancement already for the first emission whereas the BFKL 
enhancement only comes in at higher order as depicted in Fig. |[ (ii) the factor giving the 
enhancement from the BFKL dynamics multiplies the cross-section both for incoming gluon 
and incoming quark, (in) the BFKL enhancement is in the hard cross-section whereas in 
the colour dipole model the enhancement comes from replacing the ratio of parton densities 
with a cut-off in transverse momentum for the emitted gluon. It should also be stressed that 
for the 2+1 jet cross-section the Ariadne prediction gives an even larger enhancement of 
the cross-section than one would expect from BFKL dynamics. 

In summary this paper has shown the underlying dynamical reason for the differences 
in 2+1 jet cross-sections, forward jet cross-sections and transverse energy flow between the 
DGLAP based Lepto Monte Carlo and the Ariadne Monte Carlo based on the colour 
dipole model for DIS. The dominating difference is the 'unorthodox' suppression factor 
used instead of the ratio of parton densities and not the pj_-ordering. The difference can 
essentially be parameterised by a factor x"*, where h is in the order of 0.5 — 1.0 for small 
X, which multiplies the cross-section for additional emissions in Ariadne. For the forward 
jet production and the transverse energy flow this factor x^^ resembles qualitatively what 
one gets from a BFKL calculation where h = UsCa'^^'^'^/'^ — 0.5. However for the 2+1 jet 
production one rather expects an enhancement from BFKL dynamics as a function of the 
rapidity difference Ay = ln(s/|t|) between the jets. Thus, it should be possible to use the 
2+1 jet cross-section for fixed ^ and its dependence on the rapidity difference as a probe for 
BFKL dynamics. 
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